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PREPARATION OF THIS DOCUMENT 

As telenet is becoming more and more important in African inland waters where 
PAO Projects are operating, new techniques have been used in a variety of instances 
for quick assessment of the characteristics of the surveyed populations. A quanti- 
tative acoustic survey carried out in Lake Tanganyika, in cooperation with the 
Pishery Projects of Tanzania and Burundi, aims to map the distribution and estimate 
the abundance of the pelagic fish in the lake. This document discusses the statis- 
tical aspect of the echo survey design and operation which govern the accuracy and 
precision of the abundance estimation. 
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ABSTRACT 



An acoustic survey!/ assessed the total biomass and distribution of pelagic 
fish in Lake Tanganyika (Lat. 5*S, Long* 29*E) in November 1973. The lake was post- 
stratified and by using graphic methods, total biomass of the surveyed population 
was estimated at about 2.8 million metric tons. Night sample observations were 
significantly higher than day-light observations with an estimated average ratio 
4:1. Night estimates were strongly over-dispersed and fitted the negative binomial 
distribution. Day estimates were also over-dispersed and after a proper smoothing 
of the empirical data fitted the negative binomial distribution. The arithmetic 
data were log-normalized for use in any further statistical analysis. 

The acoustic surveys have been compared with the traffic surveys and suggestions 
have been made for the development of Quality Check Surveys of Echo Surveys. 

The Mean Difference of GINI has proved to be useful for providing indications 
of the level of dispersion of pelagic fish in the lake by stratum and by day and 
night periods. Also, the dispersion pattern of pelagic fish in the lake can 
effectively be portrayed by use of oblique projection charts. 

The sources and causes of potential errors inherent in echo surveys are dis- 
cussed and their spectrum has been classified into domains by taking into account 
the nature of the errors. It is indicated that the joint effect of the various 
sources of errors (sampling and non-sampling errors) are responsible for the 
observed differences between night and day sample observations. 

An estimator based on the method of collapsed strata has proved effective for 
estimating both the size of the total biomass of pelagic fish and its precision. It 
has been proved that the level of efficiency of the line sample of the survey is a 
function of the density of fish, and that for increasing the precision of sample 
estimates the distance between tracks should be reduced in proportion to the den- 
sity of the surveyed fish population. 



The required data for this study were kindly provided by Mr. K.A. Johannesson, 
Fishery Resources Officer (Acoustic Surveys), FAO; see also: K.A. Johannesson, 
Preliminary quantitative estimates of pelagic fish stocks in Lake Tanganyika 
by use of echo integration methods, EIFAC/74/ I/ Symp . 54 
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CHAPTER 1: AN OUTLINE OF THE PROBLEM 

I . 1 Introduction 

In tht field of applied fishery research quick assessment of the absolute 
abundance and distribution of fish stocks is of a great theoretical and practical 
importance. Acoustic assessment techniques have been used in a variety of instances 
and are rapidly replacing the older classical methods, primarily because of the 
rapidity with which estimates can be made and the relatively low cost of operation 
beyond the initial outlay. 

1 . 2 The nature of the problem 

The large-scale acoustic survey at Lake Tanganyika (surface area of 10 000 
n mi 2 ) was apparently the largest of its kind in freshwater research, a fact that 
increases the importance of the obtained results and their interpretation. A cri- 
tical evaluation of the survey system of the survey would provide the information 
needed for an assessment of the level of efficiency of the design of the survey and 
the level of reliability of the sample estimates. Specifically, this paper services 
two purposes: (a) According to the design of the survey, estimates of the absolute 
biomass were calculated by using a "graphical" method. This method of estimation, 
although simple, cannot be used to calculate the precision of the sample estimates. 
The need for developing proper estimators for the survey is obvious; (b) Echo 
surveys should be classified in the category of "sophisticated" surveys. Although 
the quality of the field operations of the echo survey in Lake Tanganyika was high 
it is very desirable at this stage to check the validity of the measuring device. 
Specifically, by taking into account the peculiarities of the surveyed population, 
we detect the potential sources of errors inherent in the survey system of the survey 
and estimate their effect on the results of the survey. 

CHAPTER 2: THE DISTRIBUTION OF FLOW OF PELAGIC FISH 

By using imagination one can compare the flow of fish at any point in a lake 
per unit of time with traffic flow. The specimen-mileage corresponds to vehicle- 
mileage and expresses the number of miles swum by a fish in a water-road (three- 
dimensional space) or set of roads during a specified period. In order to simplify 
my analytical work I would like to supplement my terminology describing the distri- 
bution of flow of pelagic fish: 

(a) Total survey area (A): Is the total area in a lake covered by 
an echo survey. 

(b) If one can imagine that, at time t Q , all the water-roads laid 
down end to end, the corresponding total area (surface area) 
would be called the total fish area!/ (A f ). 

(c) The general "packing"^./ of fish within a water-road or set of 
water-roads can be expressed by layers of fish or shoals of 
fish. Further, a scaling system should be established 
expressing qualitatively the level of packing of fish, e.g. 
very dense, dense, etc. A quantitative assessment of the 
level of packing of fish is expressed by the quantity of fish 
per area unit (number of fish/area unit or weight of fish/area 
unit) . 



jy Although a volume-unit is more realistic than an area-unit in this context, 
to simplify our analysis one can think of fish in a three-dimensional sense 
in a lake being projected on a two-dimensional surface 

2/ From an aggregation point of view 
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It seems to sit that, in any attempt for estimating the dispersion pattern of 
a fiah population at time t o , tha paculiaritiaa of tha survey population ahould be 
taken into account. One can aaa apacific diaparaion pattarna rather than a general 
diaperaion pattern of the fiah population. For example, the general packing of one 
tub-population of fiah, the quantitative level of packing of fish, and the apatial 
diatribution of fiah in a aet of roada might be different from that of a aecond sub- 
population of fiah occupying another aet of roada. Innate behaviour pattern! of 
individuala within a ahoal and other factors might be responsible for the statis- 
tical significance between the organic atructure of the sub-populations. 

As one can see, our imaginary water-road system is not fixed over time. Also, 
on a static basis, it is impossible to use simple geographical rules in order to 
define a water-road system in a lake. If this information was available at the 
process of designing an echo survey, the efficiency of the sampling design of the 
survey would be very high. 

It haa been argued that an objective assessment of the distribution of flow of 
pelagic fish is impossible prior to the design of an echo survey. However, indica- 
tions of the distribution of flow of pelagic fish at time t Q , can be obtained after 
the survey has been conducted and by using the sample observations of a properly 
designed echo survey. 

CHAPTER 3: SAMPLE SURVEY PROCEDURE (LAKE TANGANYIKA) 
3 . 1 Composition of pelagic stocki-/ 

In Lake Tanganyika all pelagic biomass is composed of six species. The prey 
include the two small clupeids Limnothrissa miodon and Stolothrissa tanganyika. The 
four predators comprising three very similar endemic species of the genus Lates (JL. 
angustiphrons I*, mariae and L. microlepis) , and another endemic predator of smalTer 
size and presumably shorter life history, Luciolates stappersii . Catch data indi- 
cate that more than half of the catch is composed of clupeids. In the survey 
reported here no attempt was made to differentiate biomass estimates by species. 

3.2 The sampling method 

From a statistical point of view the sampling method of the echo survey at Lake 
Tanganyika can be described as a Stratified Random Line Sample (SRLS). Past ex- 
perience has proved that this method of sampling is useful when an area sampling 
frame is used for the selection of the sample of a survey!./. In our case, pre- 
determined survey tracks were designed as a parallel grid with a total of 48 tran- 
sects across practically the whole length of the lake at 6 n mi intervals for the 
northern half of the lake, and at 10 n mi intervals for the southern half. The 
survey track in Burundi/Zaire waters was run twice giving two estimates for that 
part of the lake which were separated by a 14-day intervall/. 



J_/ See also: Henderson, Ryder and Kudhongania (1973): Assessing Fisheries 
Potentials of Lakes and Reservoirs, FAO-FI : FMD/73/S-32 

_2/ See also: The Design of Fisheries Statistical Surveys, Inland Waters, FAO- 
FIPS/T.133, by G.P. Bazigos, p. 80 

3/ The field work of the survey was carried out during a 30-day period between 
20 October abd 24 November 1973, of which 12 days were spent working out 
methods to catch and preserve live fish, in order to calibrate the equipment 
experimentally against pelagic fish species 
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3 . 3 Field operational/ 

For the field operations of the survey the research vessel LADY ALICE II of the 
Burundi Fishery Project (LOA 39 ft) wee used. A Simrad scientific sounding system, 
consisting of en EK 120 (120 KHZ) sounder coupled with complementary electronic 
instrumentation, wes temporarily installed in the wheelhouse. The instruments 
included: echo integrator QM-MK II, Hewlett Packard storage oscilloscope, electronic 
counter, signal generator, AC amplifier-voltmeter, and signal attenuator* A cali- 
brated test hydrophone was also included. 

For survey purposes the vessel speed was adjusted to 6 knots and continuous 
manual data logging was carried out every 20 minutes, thus giving 2 n mi as the 
elemental sampling distance unit (ESDU) . 

Generally, the two integrator channels (A and B) , with 10 decibel difference 
in gain setting, were adjusted to cover the whole depth interval of the fish distri- 
bution. This procedure prevented vertical stratification of the echo integration, 
but it made it possible to cope with the wide dynamic range of echo intensities 
encountered . 



3.4 



The sample of the survey 



In order to simplify the field operations of the survey and provide estimates 
on a regional basis, the lake was divided into five areas, here called strata. In 
the Table below (Table 3.4.1) the established strata are given. Also, information 
is provided on the number of sample ESDUs and area of individual stratum (n mi 2 ) of 
the surveyed population. 



Table 3.4. 



The sample of the survey 



Strata 


No. of 
sample ESDUs 


Total survey area 


Name 


Location 


2 mi 


Z 


n mi 2 


Z 


TOTAL: 

1 . Burundi/ 
Cap Banza 


(Lat.) 

03*20'- 
04*30'S 


600 


100.00 


6 458 


100.00 


104 


17.33 


698 


10.81 


2. Area near 
Kigoma 


04*30'- 
05*40'S 


175 


29. 17 


1 410 


21.83 


3. Area near 
Lagosa 


05*40'- 
06'30'S 


104 


17.33 


1 121 


17.36 


4. Area near 
Kipi li/Karema 


06*30'- 
08*00'S 


160 


26.67 


2 446 


37.88 


5. Zambia 


08*00'- 
08*45'S 


57 


9.50 


783 


12. 12 



From the tabulated deta (Table 3.4.1) one can see that the allocation of the 
total number of sample ESDUs to the strata is more or less proportional to the area 
of the strata. 



I/ 



See footnote on page iv 
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3.5 Tht procedure for estimating total biomass 

According to the sampling design of the survey the estimated total fish biomass 
(Wg, weight) can ba ragardad aa tha product of two magnitudes!/: 



(I) 



5 > -Jl Vh(B) 

where, . 

A. : Area of the survey population (n mi 2 ) in the h 
stratum 

Pu/g\ 5 Estimated average fish biomaas density (t) per 
' area sampling unit (ASU - I n mi 2 ) in the h th 



tratum 



(2) 



Further, 3^(6) * tn< product of two magnitudes each of which is 
estimated independently. 



where, 



n. : Estimated average integrator readings (mm) per ESDU 

n * . r h 



in tha h 



stratum 



C : Estimated general converting factor. On the basis 
of C the integrator readings per elemental sampling 
distance unit are converted into metric tons per ASU 
(t/n mi 2 ). 

For the estimation of R the sample data of the acoustic survey were used. An 
estimate of C was obtained through an ad hoc simulation survey, i.e. calibration of 
the acoustic equipment experimentally against a standard density of live fish (see 
footnote on page iv ). 

From the foregoing discussion it is obvious that the level of reliability of W ft 
is a function of the level of reliability of ftt which in turn is a function of the 
levels of reliability of S and n. 

CHAPTER 4: A QUICK VALIDATION OF THE RESULTS OF ECHO SURVEYS 
4 . I Introduction 

In this chapter, the need for development of quality check surveys for a quick 
validation of the results of main echo surveys is discussed. 



4.2 Quality Check Surveys of Echo Surveys (QCS-ES) 

It has been discussed in a previous chapter (Chapter 2) that by using imagina- 
tion the flow of fish at any point in a lake per unit of time can be compared with 
traffic flow. In such a case, one can argue that sampling methods can be developed 
to provide reliable estimates both of the volume of fish traffic on a static basis 
and the trend of this traffic on a dynamic basis. Further, tha results of a pro- 
perly designed fish traffic survey, baaed on small samples and running concurrently 
with an echo survey, can ba used to check the level of reliability of tha results 
of tha main echo survey. For our purposes we shall call tha specific fish traffic 
survey Quality Check Survey of tha Echo Survey (QCS-ES). 

jy For integration purposes, in this paper we can use the same notation as 
employed by K.A. Johannesson; see footnote on page iv 
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It should be noted that the object of the echo survey at Lake Tanganyika was 
to provide quick estimates of the area distribution and absolute abundance of the 
pelagic fish biomass in the lake. For this purpose the sampling method of the 
survey was based on Stratified Random Line Sample (see Chapter 3). By using 
imagination one can see that the sample tracks of the survey are check points 
(lines) in the water-roads, and the sample observations indicate the volume of fish 
passing through the check points at instantaneous periods of time. As one can 
observe, the echo survey does not collect information on the volume of fish traffic. 

The justification of my theoretical proposals for using specific fish traffic 
surveys for a quick validation of the results of echo surveys can be explained in 
simple language as follows. Assuming that by using certain control characteristics 
for stratification the lake is divided into a number of vertical hydrographic zones 
and that each zone represents a set of water-roads. Assuming also that a syste- 
matic sample of check points is selected within each zone, then by using a proper 
device, information can be collected at each sample check point on the fish traffic 
for specified periods of time. The length of the "reference period" at checking 
points should be a function of the dispersion pattern and the level of movement of 
pelagic fish. As A parenthesis, I would like to say here that, in my opinion, one 
of the sources of errors responsible for the observed differences between night and 
day measurements could be attributed to the experimental conditions of the survey 
as far as the "reference period" is concerned!/. Specifically, in the echo survey 
the sample observations were based on instantaneous reference periods of time (t ) 
for day and night measurements. By doing so we introduced the assumption of 
homogeneity over time as far as dispersion pattern and level of movements of pelagic 
fish are concerned. This, however, is not the case in Lake Tanganyika, 

In the Figure below a grid of the line sample of the survey is portrayed, 
consisting of two parallel tracks (AB, CD) and the vertical line BC. 



From the foregoing discussion one can note that only the parallel lines AB, CD 
should be considered as check points in our water-road system. Measurements along 
the line BC should be excluded because they provide information on the following two 
heterogeneous magnitudes: (a) vertical distribution and abundance of pelagic fish 
moving in the same direction with the research vessel; and (b) measurements of the 
volume of fish crossing the line on an instantaneous reference period basis. 

If there was regularity in the dispersion pattern and the level of movement of 
pelagic fish, one could expect that the sample observations at AB should not differ 
significantly from those at CD, and that the t e should be considered as the proper 
reference period for the survey system of the survey. In such a case the need 
of a QCS-ES can hardly be justified. One can argue that differences in the values 
of average integrator readings per elemental sampling distance unit between the 
parallel lines (total or stratified), give an indication of the level of irregularity 
in the dispersion pattern and level of movements of pelagic fishi/. 

y See Chapter 6 
27 See Chapter 5 
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QCS-ESs are of value when the Assumption of regularity as described above is 
not valid. In such a case QCSs are needed to detect the potential sources of error 
and to estimate their effect on the results of the main echo surveys. 

External estimates, i.e. estimates obtained through classical methods of esti- 
mating standing crops of fish, can also be used in order to provide indications of 
the level of validity of the results of an echo survey. 

CHAPTER 5: ESTIMATED DISPERSION PATTERN OP PELAGIC PISH IN LAKE TANGANYIKA 

1. I Introduction 

In this Chapter, by using the sample data of the echo survey, indications are 

provided regarding the dispersion pattern of pelagic fish in the lake* To permit a 

simple presentation of the dispersion pattern, elementary statistical methods have 
been usedi-'. 

5.2 Estimated fish area (A f ) 

Are pelagic fish everywhere in the lake (survey area)!./? The answer to the 
question is no. The estimated confidence limits for the proportion of ESDUs with 
fish (p^) are for the lake as a whole 69 percent lower limit, and 80 percent upper 
limit (a 5 percent). It is worth noting that, with the exception of Str.l, there 
are no significant differences in the estimated values of p^ between- the individual 
stratum and for day and night measurements within strata. The Table below (Table 
5.2.1) gives the estimated relative distribution of the survey area into fish area 
(p^), and the remaining area (q^ 1~P{) by stratum and by day and night periods of 
time. Also, Table 5.2.2 gives the estimated confidence intervals for P (a 0.05). 



Table 5.2.1 Estimated percentage distribution of the survey 
area into fish area (p^) and the remaining area 
(q^) by stratum and by day and night periods of 
time 



Day time/proportions 


Space 


Propor- 
tion 


Total 


Str. 1 


Str. 2 


Str. 3 


Str. 4 


Str. 5 


Total: 

Pish 
area 

Remaining 
area 


p i 
q i 


1 .0000 


1 .0000 


1 .0000 


1 .0000 


1 .0000 


1.0000 


0.7457 
0.2543 


0.7500 
0.2500 


0.8072 
0. 1928 


0.6591 
0.3409 


0.7812 
0.2188 


0.3333 
0.6667 


Night time/proportions 


Total: 

Pish 
area 

Remaining 
area 


Pi 
q i 


1.0000 


1 .0000 


1 .0000 


1.0000 


1.0000 


1.0000 


0.6973 
0.3027 


0.9200 
0.0800 


0.6167 
0.3833 


0.8000 
0.2000 


0.7500 
0.2500 


0.4167 
0.5833 



jy See also Allendix I, Mathematical Model for Acoustic Survey in Lake Tanganyika 
2/ Por our analysis the sample data of the horisontal tracks only were used 
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Table 5.2.2 Estimated confidence liniti for P, 
(P - 95Z) 



Day tim 


Total 


Str.l 


Str.2 


Str.3 


Str.4 


Str.5 


Win, 


Max. 


Min. 


Max. 


kin. 


MAX. 


Min. 


MAX. 


'Sin. 


MAX. 


Min. 


MAX. 


0.6895 


0.8019 


0.5976* 


0.8024* 


0.7218 


0.8926 


0.5174 


0.8008 


0.6791 


0.8833 





0.7176 


Might tU* 


0.6307 


0.7639 


0.8030* 


1.0000* 


0.4926 


0.7408 


0.6544 


0.9456 


0.6206 


0.8794 


0.2136 


0.6198 



* The samples have been selected from two different populations 
at far as the surveyed characteristic is concerned 



Figure 5.2.1 provides a graphical explanation of the relative composition of 
the survey area. As one can see, an estimate of the Absolute size of the fish Area 
(Af) can be calculated by applying the estimated proportions (p*^) on the total survey 
area which is a known magnitude (see Chapter 7). 



5 . 3 Estimated level of dispersion 

In our case the Mean Difference of GINll/ (d f ) can be used in order to provide 
indications of the level of dispersion of pelagic fish in the lake. Specifically, 
GINI's method can be employed to provide both estimates of the absolute differences 
between one unit and another of the survey variate and a measure of the level of 
dispersion of the respective empirical distributions (average of the absolute dif- 
ferences between one unit and another). 

By using the sample data I prepared the respective matrices indicating the 
absolute differences between the individual values of ESDUs (integrator's readings, 
mm * see Appendix II). Specifically, calculations were made on a stratum basis and 
separately for day and night periods of time (0600-2000, 2000-0600). 

In the Table below (Table 5.3.1) the estimated values of d l are given by stratum 
and by day and night period si/. 



Table 5.3.1 Estimated values of dj and average integrator readings per 
ESDU (y - M) by stratum and by day and night periods 



Space 


$! values (mm) 


(2) 


y(-ff) 


HI 

(6) 


Day 
(1) 


Night 
(2) 


(3) 


Day 
(4) 


Night 
(5) 


Str. 1 


27.65 


61.86 


2.24 


18.13 


48.83 


2.69 


Str.2 


39.61 


197.33 


4.98 


25.54 


1 16. 10 


4.55 


Str.3 


15.75 


255.36 


16.21 


10.61 


133.53 


12.59 


Str. 4 


150.42 


487.41 


3.24 


91.98 


284.20 


3.09 



\l 



See: G.P. Bacigos - Applied Fishery Statistics, FAO-FIPS/T. 1 35 , P. 33-36 



2V The sample observations of Stratum 5 were insufficient for this kind of 
analysis 
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Figure 5*3.1 Trends in space of the magnitude! 
(confidence interval P 95Z) and 
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The tabulated data (Table 5.3.1) indicate that the estimated level of diaper- 
aion of pelagic fish expressed by the valuee of 3 t ia higher during the night period 
than during the day period, and that there ia a positive relationship between the 
abundance of fish (y) and level of diapersion of fish (dj). 

Another point of interest in this kind of analysis is the high variability in 
space of i t , column 6 in Table 5*3.1 (& 2 estimated ratio of the abundance of fish 
between night and day periods). Other things being equal, based on two samples 
selected from the same population, one would expect that if the values of p^ do not 
differ significantly in space (between strata), then the eatimated values of 2 
should follow a regular pattern. This, however, is not the case in strata 2,3, and 
4, in which the hypothesis, Ho : P f - P f P % is valid. A graphical presentation 
of the estimated trenda in space of the magnitudes p* and . ia given in Figure 
5.3.1. 

5 . 4 A graphical presentation of the dispersion pattern 

It is a common practice in the field of applied fishery statistics to use charts 
for presenting and explaining numerical data. In our case, charta in oblique pro- 
jection are used to portray the dispersion pattern of pelagic fish in Lake Tanganyika. 

A first assessment of the sample data of the echo survey indicates that, very 
roughly, we can establiah three aets of water-roads in the lake: 

Water-road I - (Lj): It covers the east littoral cone of the lake 

Water-road 2 - (L 2 ): Open watera 

Water-road 3 - (L,): It covers the west littoral zone of the lake 

Figures 5.4.1 to 5.4.10, prepared for day and night periods on a stratum basis 
and for the lake as a whole, are baaed on the aample average integrator readings per 
ESDU (y). The diagrams allow an optical comparison of the magnitude y vertically, 
between tracks within strata and between strata for the lake as a whole, and hori- 
zontally between water-roads. The values of the bars of L l can be simply read on 
the y-axis, whereaa the values of the bars of L 2 and L, can be read on the y-axis 
only by taking into account the correspondingly higher drawn axes of L and L . 
The numbers in the charts are the values of the estimated specific indices expressing 
relative changes of y between the establiahed water-roada. 

A visual assessment of the prepared chart-units leads to the following 
conclusions: 

(a) Night measurements: There is consistency in the portrayed 
distribution of fish in space and L 2 can be considered as 
the modal area of the distribution. 

(b) Day measurements: The consistency of the portrayed dis- 
tribution of fish in space is a function of the recorded 
abundance of fish. Low abundance of fish fails to por- 
tray any clear dispersion pattern of fish (Str.3), 
whereas high abundance portrays a distribution similar 
to that of the night period!/ (Str.4). 



See alao: Chapman, D. (1974) - Distribution of biomass of pelagic fish in 
Lake Tanganyika, FAO-URT/7 1 /O I 2/7 . Working Paper No. 11 
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CHAPTER 6: TYPE AND SOURCES OF ERRORS XM ECIO SURVEYS!/ 

1. I Introduction 

In this chapter we diacuas tha varioua types of errors (sampling arrora and 
non-sampling arrora) inharant in tha aurvey system of an acho survey and thair affact 
on tha results of tha survey. It ia believed that tha direction and aiaa of tha 
various typaa of arrora of tha aurvay (ona occaaion) ia a function of tha dial 
behavioural pattern of the surveyed population* 

6.2 Sources and cauaea of errora 

From a sampling point of view, detailed knowledge of the possible aourcea and 
causes of errors of a large-scale survey will help to minimite their operation through 
the adoption of appropriate survey procedures and techniquea. It ia safe to assume 
that in echo surveys certain factors have been left free to operate and form errora 
in the results of the surveys. Looking at the data of the main echo aurvey in Lake 
Tanganyika from the point of view of diatortion cauaed by the error factors* I pre- 
pared the block diagram given below (Figure 6.2.1). I have divided the spectrum of 
potential errors into two domains, i.e. (1) External errora (of statiatical sampling 
nature); and (2) Internal errors (non-sampling errors and technical errora). 
Specifically, external errora are due to the level of efficiency of the line sample 
of the survey which is a function of tha type of apatial distribution of the aurveyed 
population (see also section 6.3). Internal errora are attributed to the beheviour 
pattern of the surveyed population and measuring procedure of the survey. 

6.2.1 Internal errors in echo surveys 

The main sources of non-sampling errors and technical errors in echo surveys 
have been grouped into four categories by taking into account the nature of the source 
of errors (Figure 6.2.1), i.e: 

2.1 Coverage errors of type-I 

2.2 Coverage errors of type-2 / Avoidance errors 

2.3 Measurement errors 

2.4 Other kind of errors 

Further, within each of the above main sourcea of errora, I recorded the various 
factors which operate and form the respective errors (Figure 6.2.1). Specifically: 

2.1 The coverage errors of type-1 originate primarily from 
limitations regarding the minimum depth et which the 
acoustic beam produces quantifyable fiah echoea. This 
depth usually ranges from 5-7 m and ia determined by 
two factors: (i) The depth location of the transducer 
which, for hull mounted transducers, may vary between 
2 and 4 m, depending on veaaal type; and (ii) The 
feature of the T.V.G. (time varied gain) function of 
the receiving ayatam which precludea an efficient 
aignal receiption for an initial time period corres- 
ponding to a range of 3 m immediately below the 
transducer. 

Secondly, there are coverage errors (usually not very 
serious when dealing with pelagic speciea) which depend 
on the beam shape and pulae length employed, plus the 
topography of tha aaa bed. 



I would like to acknowledge that I received very uaeful comments from K.A. 
Johannesson (FAO-Fishery Resources Officer) during the preparation of thia 
chapter 
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2.2 The magnitude of the coverage errors of type-2/ 
avoidance errors it apparently a function of a 
number of factors, aucn as: 

(i) The intensity of the noise field radiated 
by a moving vessel 

(ii) The vulnerability of the fish species and 
hence reaction (fright?) to this noise, 
which can be conceived as a disturbance- 
cone of a high level near the surface but 
diminishing with increased depth 

(iii) The degree to which this zone penetrates 
the natural environment of the species 
concerned, e.g. those which typically 
ascend and form near-surface shoals 
during day time, or species that always 
remain relatively near the surface simply 
because of shallow water depths 

2.3 The measurement errors are attributed to the inaccuracy 
of the electroacoustic measuring device and their overall 
magnitude largely depends on the accuracy with which 
system calibrations can be performed. 

2.4 Other kind of errors are those related to: 

(i) The fact that the vessel is pitching, rolling 
and tossing, etc., and the transducer (when 
stabilized) has to move with it 

(ii) The operators skill to optimize control settings 

6. 3 Differences between day and night observations 

Directly connected with the effect of the various types of errors described 
above, are the observed differences between day and night measurements,!/ . In spite 
of the standardization of the field operations of the survey, the magnitude of 
external errors is higher during the day time rather than during the night. Imagine 
someone at a cross-road, without traffic lights, and he crosses the road straight 
away (instantaneous reference period) without taking any precaution. From a mathe- 
matical point of view, the probability (P) of the person crossing the road being hit 
by a car is a function of the level of flow of cars. If there is no flow of cars, 
the value of P is equal to zero. If there is a continuous flow of cars the value 
of P is equal to 1, 0<P<1. Furthermore, if the traffic flow per unit of time is 
given the probability of a person crossing the road and being hit by a car is a 
function of the "packing 1 * pattern of cars. This in turn means that with an instan- 
taneous reference period in the survey system of our survey, the level of efficiency 
of the line sample is a function of the packing pattern of moving fish. The effi- 
ciency of the sample is higher when the packing pattern of fish is expressed by 
layers of fish (night time observations) rather than when it is expressed by shoals 
of fish (day time observations). External errors lead to an under-estimation of 
the surveyed magnitude in day time. 

It has been observed that during the survey period fish biomass tended to con- 
centrate closer to the surface during the day rather than during the night. This in 
turn means that the magnitude of the internal errors of type 2.1 and 2,2, i.e. 
coverage errors of type-1 and coverage errors of type-2 and avoidance errors are 
bigger in day time rather than in night time. These errors lead to an under- 
estimation of the surveyed magnitude in day time. 



\J The estimated average overall ratio between night and day measurements was 4:1 
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Finally, it hat been indicated that tha maaauring accuracy of tha device ia 
battar whan wa ara dealing with layera of fiah (night time maaauramanta) rathar than 
with ahoala of fiah (day tima measurements) . 

From tha foragoing discussion ona can aaa that tha joint affact of tha various 
aourcas of arror (intarnal and aztarnal) ara raaponaibla for tha obaarvad diffarancas 
batvaan day tima and night tima maaauramanta. 

To gat an idaa of tha joint affact of tha varioua aourcaa of arrora on tha 
accuracy of tha sampla aatimataa, I praparad Figura 6.3*1. In tha diagram!/ thara 
ara graphical praaantationa of tha confidanca intarvals (P 95Z) for ? h (population 
avaraga intagrator raadinga par ESDU) for day and night and within tha aatablishad 
strata. A viaual assaaamant of tha respective charta indicataa that: 

Str.3 Thara ara highly aignificant diffarancas between 
day and night obaarvations 

Str.2 Thara ara significant diffarancas batween day and 
night obaarvationa 

Str.1,4 Tha astablishad confidanca intarvals are overlapping 

Tha observed inconsistency in tha level of atatistical aignificance between 
day and night estimates of the surveyed magnitude portrays the unknown weighting 
system on the basis of which the various sources of errors affect the sample 
observations . 



\J For calculations the sample data of the horizontal tracks only were taken 
into account 
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CHAPTER 7: ESTIMATION OF POPULATION VALUES 
7. 1 Introduction 

In this chapter I determine the estimators which can be used to calculate the 
surveyed population values (eee Figure 7.2*1) and their sampling errors!'. Also, 
the relative efficiency of the different estimation procedures is discussed. 

7.2 Method of collapsed strata 

One method which can be used to estimate the population values is the method of 
collapsed strata* For the application of the method we should extend the stratifi- 
cation of the lake to the point where the sample contains only one track in each sub- 
stratum. For our purposes, and in order to simplify the analysis of the available 
sample observations 9 a sub-stratum is formed by considering the water up to a miles 
on either side of a given sample track/, i.e. a equals three miles for the northern 
part of the lake, and five miles for the southern part. In this method, in which a 
track has been selected from each sub-stratum, the sampling variance is estimated by 
pairing the sub-strata to form collapsed strata* Since the precision of the method 
increases if the pairing can be so arranged (before the collection of data) that the 
strata forming the pair are about equal in size (total or average of the survey 
variate), I formed collapsed strata by pairing adjacent sub-strata within the 
existing strata. 



Notation!/: 



i 
k 

y hi 
n hi 



stands for a given sub-stratum 
stands for a given collapsed stratum 

total integrator reading (mm) of the 
sample track in the i ch sub-stratum, 
within the h th stratum 

number of ESDUs in the sample track 

: Estimated total biomass (t) 
in the i sub-stratum 

: Estimated total biomass (t) 
in the h stratum 

: Estimated total biomass (t) 
for the lake as a whole 



J_/ Estimators have been developed under the following two assumptions: (a) in the 
field operations of the survey, the same shoals/layers have not been covered 
more than once over the successive sample transects of the survey (see also 
Chapter 4); and (b) the calibration method of the survey was free of errors. 
The statistical efficiency of the calibration method of echo surveys will be 
discussed in another paper* Also, in the same paper, I will determine the 
estimator of the overall variance of the sample estimates and its analysis 
into components. For estimation purposes the night sample observations of 
the horizontal tracks of the survey were used. 

2/ Under ideal conditions the lake should be stratified in advance into a number 
of limnological zones (horizontally) and one sample track should be selected 
randomly within the established sub-strata (limnological zones). 

3/ See also section 3.5 
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i Estimated varianca of tha estimated 
total biomass in tha h t stratum 

: Estimated varianca of tha eatimated 
total biomaas for tha laka at a whole 
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*<Ak> - 


<0 5 h. 


_ S \i 
8*h2' 


t Katiaatad variance of the aatimatad 
total bioaaaa in tha k ch eollapaad 
tratuml/ 



By applying tha estimator! aatahliahad above on tha sample data estimates were 

calculated of tha survey characteristics, Tha Tabla below (Tabla 7.2.1) shows both 
tha estimates of tha total biomass of palagic fiah by stratum, major atratum and for 
tha lake as a whole, and of tha aampling variance, sampling error and relative samp- 
ling error of the estimated magnitudaa. Alao, in columns 7 and 8 of tha Table, the 
confidence limits (P 95X) of tha paramatara undar estimation are given!/. Major 
strata were formed by grouping together tha aatabliahad strata. Specifically, 
Major Stratum I consists of Strata 1,2,3 (northern part of tha laka with sample 
transects at 6 n mi intervals) and Major Stratum II conaists of Strata 4,5 (southern 
part of the laka with sample transacts at 10 n mi intarvala) . 

Tabla 7.2.1 Point astimataa and estimated confidence intervals 
for tha total biomass of palagic fiah in tha laka 
(P 95X) by atratum, major atratum, and for the 
laka as a whole (in metric tons) 



Maj or 
Str. Str. 


Point 
astimataa 


Sampling 
variances 


Sampling 
arrora 


Relative 
sampling 
arrora 


Conf. intervals 
(95X) 


Lower 
limit 


Upper 
limit 


(1) (2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


TOTAL: 


3051953 


438664796780 


662317 


21.70 


1612698 


4491 168 


I 
1 


997717 


22145267449 


148812 


14.92 


524198 


1471236 


88686 


104572816 


10226 


11.53 


56147 


121225 


2 


474930 


17039000997 


130533 


27.48 


a 


* 


3 


434101 


5001693626 


70722 


16.29 


* 


* 


II 
4 


2054236 


416519529331 


645383 


31.42 


* 


* 


2015940 


416256862482 


645179 


32.00 


* 


* 


5 


38296 


262666849 


16207 


42.32 


* 


* 



The sampling distribution of tha estimate is highly dispersed 



The quantity v (gW h| ) 
The overstatement 



overstates tha true varianca of gW... 
depending upon tha extant to which 
strata forming tha same pair differ with raapact to their totals 
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One should ramambar that tha siie of tha total biomass is not fixed over time 
and that tha tabulated data (Tabla 7.2.1) gives tha sise of tha estimated 
magnitudes at a givan period of time, i.e. November 1973 
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An assessment of the results of the method of collapsed etraTr*<Table 7.2.1) 
leads to tht following conclusions: 

(a) Tht strata variance is a function of tha laval of abundance of 
pelagic fish and proportional allocation of the sample tracks 
between strata is not a beneficial method. In our case the 
optimum allocation must be used. 

(b) The sample sise in Major Stratum II (sample transects at 10 
n mi intervals) should be considered as a small sample for 
the survey and this is reflected in the estimated relative 
sampling error. 

(c) The level of precision achieved in Major Stratum II has 
affected the level of reliability of the sample estimates 
of total biomass of pelagic fish for the lake as a whole. 

(d) The level of precision of the estimated magnitude for Major 
Stratum I should be considered as good, 

(e) The precision of the method of collapsed strata would be 
highly improved if the sample of the survey in the strata 
with a high abundance of fish (specifically Str. 2 and 4) 
was based on transects at 3 n mi intervals and in the 
remaining ones at 6 n mi intervals. 

The Table below (Table 7.2.2) gives the point estimates based on the methods 
of collapsed strata and of the "graphical method 1 * used in the turveyl/. 



Table 7.2.2 Point estimates of the total biomass of pelagic 
fish in Lake Tanganyika based on various 
estimation procedures (in metric tons) 



Major 
Str. Str. 


Point estimates 


Point estimates (Z) 


Graphical 
method 


Method of 
collapsed 
strata 


Graphical method 


Method of 
collapsed strata 


TOTAL: 


2600000 


3051953 


100.00 




100.00 




I 

1 
2 
3 


820000 


997717 


31 .54 


100.00 


32.69 


100.00 


120000 
360000 
340000 


88686 
474930 
434101 




14.63 
43.90 
41 .47 




8.89 
47.60 
43.51 


II 

4 
5 


1780000 


2054236 


68.46 


100.00 


67.31 


100.00 


1700000 
80000 


2015940 
38296 




95.51 
4.49 




98.13 
1.87 



\J For estimation purposes in the main survey, the day time sample data were 
"" multiplied by 4 (overall correcting factor) before they matched with the 
night time sample data 
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Figure 7.2.1 Relative Abundance of pelagic fish stocks in 
Lake Tanganyika 
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In tha above Table the point estimates baaed on tha two estimation procedures 
do not show any statistically significant differences at P - 95Z. 

In Figure 7.2*1 tha qualitative raaults of tha echo aurvey basad on tha 
graphical method ara portrayed* 

7.3 A second estimator for total biomass of pelagic fish 

Another method to estimate tha total biomaes of pelagic fish in the lake can 
be developed by considering the sampla tracks (horisontal tracks) within the estab- 
lished strata as units of random samples selected from the survey peculation. Also, 
since we ara dealing with sampla tracks of various sites, as far as the number of 
ESDUs is concerned, the following estimator can be used for tha calculation of the 
surveyed population total: 

Notation: j : Tha suffix j stands for a given track 

h : The auffix h stands for a given stratum 

y. . : Total integrator reading (mm) of the j th 
J track in stratum h 

n. : Number of sample tracks in the h stratum 
n 

: Estimated total biomass in the h fc 
stratum based on the j th track and 
n. . number of ESDUs in the j tn track 

fi W, -^ g^hm : Efti " d total average biomass in 
p h JP IJ ' the h tn stratum 

W - I : Estimated total biomass for the lake 
" as a whole 



v * h/ vv^ 

: Estimated variance of the estimated 
total biomass in the h stratum 

v( B^ " v *B^h^ : Eftimated variance of the estimated 

total biomass for the lake as a whole 

An assessment of the calculated precision of the estimated magnitude indicates 
that the method is not the proper one for our purposes. The estimated relative 
sampling error for the lake as a whole was of tha level of about 50 percent. 

7.4 A third estimator for total biomass of pelagic fish 

An assessment of the dispersion pattern of pelagic fish in the lake (Chapter 5) 
indicates that the precision of the estimates basad on unstratifiad sample (section 
7.3) can be highly improved if the method of stratified sampling based on inter- 
penetrating sub-samples is used. In such a case the lake should be divided into 
not less than 30 strata by horisontal partition lines and into 3 main sets of water- 
roads, dividing the lake lengthwise. Two sampla tracks should be selected within 
the established cells (systematic selection, in order to simplify the field 
operations). It is expected that if tha proper control characteristics for 
stratification should be employed, this method would provide much more precise 
estimates than those of the method described in section 7.3. 
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7.5 The relative efficiency of precision 

In sampling techniques, the relative efficiency of precision of two different 
methods of sampling based on the same type of sampling unit may be defined as the 
reciprocal of the ratio of the sampling variances of the estimates given by the two 
methods when the same number of units are taken. In our case estimates of the pre- 
cision of the survey characteristic were calculated by using the method of collapsed 
strata (section 7.2) and the method of unstratified sampling (section 7.3). An 
estimate of the relative efficiency of the two methods for the lake as a whole is 
given by: 

. A. 
cv 1 



RE (pr) - i^ - T^TT - 5.31 



where, cv ( ft W) : la the coefficient of variation of the estimated 
2 magnitude based on method 2 (section 7.3) 

cv ( fl W) : Is the coefficient of variation of the estimated 



fi 



magnitude based on method I (section 7.2) 



An assessment of the estimated value of RE, r \ indicates that the method of 
unstratified sample provides five times the maximum variance of the estimated 
characteristic (total biomass of pelagic fish in the lake) as the method of collapsed 
strata. 
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A MATHEMATICAL MODEL FOR ACOUSTIC SURVEY 
IN LAKE TANGANYIKA 
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A MATHEMATICAL MODEL FOR ACOUSTIC SURVEY 
IN LAKE TANGANYIKA!/ 



Introduction 

The main objective of this analysis was to attempt to fit a mathematical model 
to the empirical data of the survey in order to explain the spatial distribution of 
the fish stocks, and to obtain a suitable transformation to normalize the frequency 
distribution of the data in order to perform further statistical analysis such as 
regression and analysis of variance. 

Methodology.!/ 

For this analysis Stratum 2, the largest stratum was selected. The basic 
variable considered was y integrator reading. Also, it was recognized that the 
fish move in schools at day time and in layers at night. Hence day time and night 
time data were treated separately. 



Night data 

Firstly, the night data were arranged in a frequency distribution table. The 
sample mean (y) and variance (s 2 ) were then computed. The purpose was to compare 
the magnitudes of the two statistics y and s 2 , to indicate which theoretical distri- 
bution might be a suitable model. Now for this type of population, three mathema- 
tical distributions are possible models, namely Positive Binomial (a 2 <y) , Poisson 
(a 2 - y) , and Negative Binomial (a 2 >u) , where y, a 2 are population mean and variance 
respectively; y was found to be 170.90, whilst s 2 was 141896.61. As ? and s 2 are 
unbiased estimates of the population mean (u) and variance (a 2 ) and s 2 >y, there is 
a strong indication that the negative binomial distribution is the most probable 
model. The chi-squared "goodness-of-f i t" test was applied to ascertain the ade- 
quacy of the model, but unfortunately the relevant negative binomial distribution 
tables were not available. Assuming the model to be adequate, the spatial dis- 
persion of the population could be described as a contagious distribution or clumped 
distribution. The negative binomial distribution is given by p (l~q) where 

"if"' ^ * *~ p and k * 8 e8t ^- mated by i ' The individual terms are given by: 



(x) 



(k+x-m 



is the probability of x individuals in a sampling unit. In our case, 
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Prepared by S. Adjei, FAO Fellow, Department of Fisheries, who was working 
under my supervision during the period 20 June to 12 July 1974 to round out 
two-year fellowship studies in the United States, Canada and Italy 

See also: G.P. Bacigos - Applied Fishery Statistics, FIPS/T.135, p. 42-50 
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Logarithmic transformation 

The next step was to find a auitablt transformation to normalize the frequency 
distribution. Now for a negative binomial distribution with k less than two, the 
usual transformation is the logarithmic transformation. In particular the zero 
counts in the data suggested the transformation Z - log t (y+1) as the appropriate 
transformation for normality. To test the validity of this transformation the chi- 
squared "goodness-of-f it" test was applied. The transformed data were assumed to 
follow the normal distribution whose mean and variance were estimated by the corres- 
ponding sample values. The expected frequencies (f|) were obtained by considering 
the cumulative distribution tables of the normal distribution. The chi-squared 
statistic, 



was then computed to be 28.11 9 where f^ - observed frequencies, and n - number of 
frequency classes. Now, under the hypothesis of normality the above statistic 
follows a chi-squared distribution with n-3 degrees of freedom (df ) . Now n 32, 
therefore df 29. But the tabulated value for 95 percent level of X 2 (?9) " 42.56. 
Hence, it is concluded that at 95 percent confidence level the normal 
distribution fits our empirical data adequately. Hence the recommended trans- 
formation is Z - log. (y+1). The analysis was repeated for the day time data and 
the negative binomial distribution was again found to be the appropriate mathematical 
model. Again the transformation log (y+1) was found to normalize the data. Except 
for two points out of 41 which could be treated as outliers all the remaining 39 
points in the frequency distribution table satisfied the chi-squared criterion. 
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STRATUM 2 - NIGHT DATA 



Mean y - Z - ^jT^ " 170.90 



2 E(fy 2 )-yfy 
Sample variance s* - ^ n -l 



687336UOO-II97496.30 
40 

141896.61 



s >y negative binomial model 



k - -4 - 0.21 

-7 



Log transformation z log 

s.d. /s 7 - s - 376,69 



Transformed data: 



log (y+1) 
Z - I = - 3.57 



n-1 
- /3.53 - 1 .88 
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ACOUSTIC SURVEY - LAKE TANGANYIKA 
STRATUM 2 - NIGHT DATA 



f 


y 


y 2 


fy 


fy 2 


Z- log^y+l) 


Z 2 f 


Zf 


1 














0.00 


0.00 


0.00 


2 


2 


4 


4 


8 


1 . 10 


2.42 


2.20 


2 


4 


16 


8 


32 


1.61 


5. 18 


3.22 




5 


25 


5 


25 


1.79 


3.20 


1.79 




6 


36 


6 


36 


1.95 


3.80 


1.95 


2 


7 


49 


14 


98 


2.08 


8.33 


4. 16 


2 


10 


too 


20 


200 


2.40 


1 1 .52 


4.80 




1 1 


121 


1 1 


121 


2.48 


6. 15 


2.48 




U 


144 


12 


144 


2.56 


6.55 


2.56 




16 


256 


16 


256 


2.83 


8.01 


2.83 




20 


400 


20 


400 


3.04 


9.24 


3.04 


2 


28 


784 


56 


1568 


3.37 


22.71 


6.74 




29 


841 


29 


841 


3.40 


11.56 


3.40 




30 


900 


30 


900 


3.43 


1 1.76 


3.43 




36 


1296 


36 


1296 


3.61 


13.03 


3.61 


2 


38 


1444 


76 


2888 


3.66 


26.79 


7.32 


2 


42 


1764 


84 


3528 


3.76 


28.14 


7.52 


2 


48 


2304 


96 


4608 


3.89 


30. 13 


7.78 




52 


2704 


52 


2704 


3.97 


15.76 


3.97 




55 


3025 


55 


3025 


4.03 


16.24 


4.03 




71 


5041 


71 


5041 


4.28 


18.32 


4.28 




72 


5184 


72 


5184 


4.29 


18.33 


4.29 




94 


8836 


94 


8836 


4.55 


20.70 


4.55 




105 


1 1025 


210 


22050 


4.66 


43.42 


9.32 




304 


92416 


304 


92416 


5.72 


32.72 


5.72 




310 


96100 


310 


96100 


5.74 


32.95 


5.74 




366 


133956 


366 


133956 


5.91 


34.93 


5.91 




450 


202500 


450 


202500 


6.1 1 


37.33 


6. 1 1 




470 


220900 


470 


220900 


6. 15 


37.82 


6. 15 




880 


774400 


880 


774400 


6.78 


45.97 


6.78 




1 170 


1368900 


1 170 


1368900 


7,06 


49.84 


7.06 




1980 


3920400 


1980 


3920400 


8.44 


71 .23 


8.44 


Tot* 


1 n - 


41 








688.08 


153. 18 



Log transformation: 



Z - log e (y+l) 
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ACOUSTIC SURVEY - LAKE TANGANYIKA 
STRATUM 2 - NIGHT DATA 



*i 


f i 


X 


F i( x) 


nF.(x) 


n(F i!r F i> 


<f t -fi> 


x 2 





1 


-1 .90 


0.03 


1.23 


1.23 


0.05 


0.04 


2 


2 


-1 .31 


0. 10 


4. 10 


2.87 


0.76 


0.26 


4 


2 


-1.04 


0. 15 


6. 15 


2.05 


0.0025 


0.00 


5 


1 


-0.95 


0. 17 


6.97 


0.82 


0.03 


0.04 


6 


1 


-0.86 


0. 19 


7.79 


0.82 


0.03 


0.04 


7 


2 


-0.79 


0.21 


8.61 


0.82 


1 .39 


1 .70 


10 


2 


-0.62 


0.27 


1.07 


2.46 


0.21 


0.09 


1 1 


1 


-0.58 


0.28 


1 .48 


0.41 


0.35 


0.85 


12 


1 


-0.54 


0.29 


1 .89 


0.41 


0.35 


0.85 


16 


1 


-0.39 


0.35 


4.35 


2.46 


2. 13 


0.87 


20 


1 


-0,28 


0.39 


5.99 


1 .64 


0.41 


0.25 


28 


2 


-0. 1063 


0.45 


18.45 


2.87 


0.76 


0.26 


29 


1 


-0.0904 


0.46 


18.86 


0.41 


0.35 


0.85 


30 


1 


-0.07 


0.47 


19.27 


0.41 


0.35 


0.85 


36 


1 


0.02 


0.51 


20.91 


1.64 


0.41 


0.25 


38 


2 


0.05 


0.52 


21.32 


0.41 


2.53 


6. 17 


42 


2 


0. 10 


0.54 


22. 14 


0.82 


1 .39 


1 .70 


48 


2 


0. 17 


0,57 


23.37 


1.23 


0.59 


0.48 


52 




0.21 


0.58 


23.78 


0.41 


0.35 


0.85 


55 




0.25 


0.60 


24.60 


0.82 


0.03 


0.04 


71 




0.3776 


0.64 


26.24 


2.05 


* 1 . 10 


0.54 


72 




0.3829 


0.65 


26.65 


0.41 


0.35 


0.85 


94 




0.5212 


0.70 


28.70 


2.05 


1 . 10 


0.54 


105 




0.58 


0.72 


29.52 


0.82 


1.39 


1 .70 


304 




. 14 


0.87 


35.67 


6. 15 


26.52 


4.31 


310 




. 1542 


0.88 


36.08 


0.41 


0.35 


0.85 


366 




.24 


0.89 


36.49 


0.41 


0.35 


0.85 


450 




.35 


0.91 


37.31 


0.82 


0.03 


0.04 


470 




.3723 


0.92 


37.72 


0.41 


0.35 


0.85 


880 




.71 


0.96 


39.36 


1 .64 


0.41 


0.25 


1 170 




.86 


0.97 


39.77 


0.41 


0.35 


0.85 


1980 




2.59 


1.00 


41 .00 


1 .23 


0.05 


0.04 



y + l x 

No. of frequency classes - 32 
No. of parameters estimated from N(p,a 2 ) 
Therefore, d.f. - 32-2-1 29 



f 2 - 28. 1 1 



Tabulated value of 



at 95 percent level - 42.56 



X (29) 

Hence normal distributions fit adequately at 95 percent 
confidence level 
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MATHEMATICAL MODEL POR ACOUSTIC SURVEY OF LAKE TANGANYIKA 

STRATUM 2 * DAY DATA 



Frequency distribution of integrator value (y) 



f 


19 


18 


5 


6 


2 


2 


1 


2 


5 


1 


1 


4 


2 


y 








2 


3 


4 


5 


6 


7 


8 


9 


11 


12 


13 



f 


i 


2 


2 


3 


3 


1 


1 


1 


2 


1 


3 


1 


1 


y 


14 


17 


18 


22 


24 


26 


30 


31 


32 


36 


42 


45 


46 



f 


2 


1 


2 




1 


1 


1 


1 


1 


1 




1 




y 


30 


54 


56 




105 


125 


133 


152 


170 


210 




360 





y 

"y? 



102 

2369 

23.23 

306091 

55031.87 



v , r(y) 



. 306091-55031.87 
101 

2485.74 



i.e. 2 >y - Negative binomial model 

k - 0.22 Required trans- 
formation y*log (y+l) 
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ACOUSTIC SURVEY - LAKE TANGANYIKA 
STRATUM 2 - DAY DATA 



*i 


f . 


X 


F.-Cx) 


nF^x.) 


*l m 
< F i + r F i> 


i-'i> f 


x 2 





19 


-1.26 


0.10 


10.20 


10.20 


77.44 


7.59* 


1 


18 


-0.82 


0.20 


20.40 


10.20 


60.84 


5.96 


2 


5 


-0.55 


0.29 


29.58 


9.18 


17.47 


1.90 


3 


6 


-0.37 


0.36 


36.72 


7.14 


1.30 


0. 18 


4 


2 


-0.23 


0.41 


41.82 


5. 10 


9.61 


1.88 


5 


2 


-0. 11 


0.46 


46.92 


5. 10 


9.61 


1.88 


6 


1 


-0.01 


0.50 


51 .00 


4.08 


9.49 


2.33 


7 


2 


0.08 


0.53 


54.06 


3.06 


1. 12 


0.37 


8 


5 


0. 15 


0.55 


56. 10 


2.04 


4.24 


2.08 


9 


1 


0.22 


0.59 


60. 18 


4.08 


9.49 


2.33 


1 1 


1 


0.34 


0.63 


64.26 


4.08 


9.49 


2.33 


12 


4 


0.39 


0.65 


66.30 


2.04 


3.84 


1.88 


13 


2 


0.44 


0.67 


68.34 


2.04 


0.00 


0.00 


14 


1 


0.48 


0.68 


69.36 


1.02 


0.00 


0.00 


17 


2 


0.60 


0.73 


74.46 


5. 10 


9.61 


1.88 


18 


2 


0.6322 


0.74 


75.48 


1 .02 


0.96 


0.94 


22 


3 


0.76 


0.78 


79.56 


4.08 


1. 17 


0.29 


24 


3 


0.86 


0.79 


80.58 


1 .02 


3.92 


3.84 


26 


1 


0.81 


0.81 


82.62 


2.04 


1 .08 


0.53 


30 


1 


0.95 


0.82 


83.64 


.02 


0.00 


0.00 


31 


1 


0.9741 


0.83 


84.66 


.02 


0.00 


0.00 


32 


2 


0.99 


0.84 


85.68 


.02 


0.96 


0.94 


36 


1 


.06 


0.86 


87.72 


2.04 


1.08 


0.53 


42 


3 


. 1612 


0.87 


88.74 


.04 


3.84 


3.69 


45 


1 


.2064 


0.88 


89.76 


.02 


0.00 


0.00 


46 


1 


.2193 


0.89 


90.78 


.02 


0.00 


0.00 


50 


2 


.27 


0.90 


91 .80 


.02 


0.96 


0.94 


54 


1 


.3225 


0.9066 


92.47 


1.02 


0.00 


0.00 


56 


2 


.3419 


0.9099 


92.81 


0.34 


2.76 


8. 1 1* 


105 


1 


.74 


0.9591 


97.83 


5.02 


16. 16 


3.22 


125 


1 


.8580 


0.9686 


98.80 


1.02 


0.00 


0.00 


133 


1 


.8967 


0.9713 


99.07 


0.27 


0.53 


1.96 


152 


1 


.98 


0.9761 


99.56 


0.49 


0.26 


0.53 


170 


1 


2.0516 


0.9798 


99,94 


0.38 


0.38 


1.00 


210 


1 


2. 19 


0.9857 


100.54 


0.60 


0. 16 


0.27 


360 


1 


2.54 


0.9945 


101 .44 


0.90 


0.01 


0.01 



z 

X 

n 



log e (yH) 



a 
102 



- o 



.55 



X 2 - 59.48 
X 2 " 43.78** 
" 2 * 43.78** 



0.05,33 



43.80 



** Estimated x 2 without the two values 
indicated by asterisk 
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APPENDIX II 



ESTIMATED DISPERSION MATRICES 
estimated Mean Difference of GINI) 
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